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ABSTRACT 

Red Supergiants (RSGs) are among the brightest stars in the local universe, mak- 
ing them ideal candidates with which to probe the properties of their host galaxies. 
However, current quantitative spectroscopic techniques require spectral resolutions of 
R^17,000, making observations of RSGs at distances greater than IMpc unfeasible. 
Here we explore the potential of quantitative spectroscopic techniques at much lower 
resolutions, R Ri2-3000. We take archival J-band spectra of a sample of RSGs in the 
Solar neighbourhood. In this spectral region the metallic lines of Fe i, Mg i. Si i and Ti i 
are prominent, while the molecular absorption features of OH, H2O, CN and CO are 
weak. We compare these data with synthetic spectra produced from the existing grid 
of model atmospheres from the MARCS project, with the aim of deriving chemical 
abundances. We find that all stars studied can be unambiguously fit by the models, 
and model parameters of log g, effective temperatures Teff, microturbulence and global 
metal content may be derived. We find that the abundances derived for the stars are 
all very close to Solar and have low dispersion, with an average of logZ=0.13±0.14. 
The values of T^s fit by the models are ^150K cooler than the stars' literature values 
for earlier spectral types when using the Levesque et al. temperature scale, though 
we find that this discrepancy may be reduced at spectral resolutions oi R = 3000 or 
higher. In any case, the temperature discrepancy has very little systematic effect on the 
derived abundances as the equivalent widths (EWs) of the metallic lines are roughly 
constant across the full temperature range of RSGs. Instead, elemental abundances 
are the dominating factor in the EWs of the diagnostic lines. Our results suggest that 
chemical abundance measurements of RSGs art possible at low- to medium-resolution, 
meaning that this technique is a viable infrared- based alternative to measuring abun- 
dance trends in external galaxies. 

Key words: stars: abundances, stars: late-type, stars: massive, supergiants, galaxies: 
abundances, galaxies: stellar content 



1 INTRODUCTION 

Red Supergiants (RSGs) are an evolved state of stars with 
initial masses 8M0<Afinit<35M0 (e.g. IMevnet fc Maeden 
[2OO0I). They have high bolometric luminosities (L ^ 
1Q*'^Lq) which peak at around l/im, meaning that they 
are tremendously bright at optical and near-IR wavelengths, 
with M,j ~ —11, rivalling the integrated light of globular 
clusters and dwarf galaxies. 

As the lifetime of massive stars is short, RSGs are neces- 
sarily young objects, with ages of ~50Myr at most. As they 
cannot have travelled very far from their place of birth, the 
chemical composition of RSGs should closely match that of 
the gas phase material in their environs, aside from a certain 



amount of CNO processing. Therefore, through quantitita- 
tive spectroscopy of RSGs it should be possible to use them 
as probes of their local chemical abundances. 

The evolutionary predecessors of RSGs, blue super- 
giants of spectral type A and B, have been very suc- 
cessfully used as tracers of chemical abundances out to 7 
Mpc distance by means of low resolu tion quantitative spec- 
tral analyis at visual wavelengths (jKudritzki et al.l I2OO8I : 
iKudritzkil I2OIC1I ) . Using individual stars in galaxies as the 
source of information about chemical composition rather 
than Hll-regions and their emission lines has clear advan- 
tages, because the l atter are subject to sub s tantial system- 
atic uncertainties |Kewlev fc EllisonI 120081 : iBresolin et al.l 
I2OO9I : IKudritzkil |201(J ). It is therefore important to find as 
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Figure 1. Spectra of RSGs in the J-band, at resolution of 2,000. Positions and identifications of the dominant spectral features have 
been indicated. 
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many reliable stellar abundances as possible useful for the 
investigation of galaxies at larger distances. The goal of this 
paper is to demonstrate that RSGs have an enormous po- 
tential in this regard. 

A difficulty in extracting chemical abundances from 
RSG spectra is that their cool atmospheres contain molec- 
ular material, and therefore their spectra are filled with lit- 
erally thousands of strong molecular spectral lines. These 
lines overlap and blend with the metallic lines. To extract 
abundance information from the spectra it is necessary to 
observe at sufficiently high resolution so as to separate the 
molcular lines and isolate the metallic features. This typ- 
ically means spectral resolutions in excess of _R > 17,000. 
Obtaining high signal-to-noise (SNR) data at this resolu- 
tion out to Mpc distances is extremely challenging, even for 
the largest telescopes. 

Abundance studies of RSGs to date have typically con- 
centrated on high-resolution observations in the _ff-band, 
where there is a rich concentration of both molecular a nd 
metallic lines fsee lCunha et al.ll2007l : iDavies et al.ll2009al lbh. 
However, if one shifts focus to the J-band, one sees that the 
molecular transitions are much reduced in both number and 
in strength, with the dominant spectral features being the 
metallic lines of Fel and the a-elements Mgl, Si I and Til. 
Hence, it should be possible to observe this spectral region 
at much lower spectral resolution while still being able to 
extract abundance information. The only drawback is that 
one loses information about stellar nucleosynthesis within 
the RSGs themselves, which is redundant information if one 
is only concerned with the abundances of the local interstel- 
lar medium. This method of deriving chemical abundances 
from J-band spectra is largely unteste d and unexplored, t he 
only such study to date being that of lOriglia et al.l (|2004l ). 

In this paper we will explore the potential of quanti- 
tative spectroscopy of RSGs in the J-band, using archival 
data and current grids of model atmospheres. In Sect. [2] we 
will describe the J-band observations of RSGs, the model 
atmospheres, and the fitting techniques we employ. The re- 
sults of our analysis are presented and discussed in Sect. O 
We conclude in Sect. \E\ 



2 DATA AND ANALYSIS 

2.1 J-band observations of RSGs 

For the majority of our analysis we use data from the IRTF 
spectral library. The spectra cover the wavelength range 0.8- 
4^m, at spectral resolution of _R = 2, 000, and at SNR>100. 
The spectra are all of bright stars in the Sol ar neighbour- 
hood. For full details of the observations, see iRavner et al.l 
l|2009l ). 

In Fig. [l] we plot these spectra, in order of increasing 
spectral type, in the wavelength range of 1.15-1.23/im. The 
spectra have been normalized by the continuum fiux. This 
continuum level is not trivial to measure in spectra with 
many features - we measured it by ranking the channel val- 
ues in the plotted range in order of decreasing flux level, 
and then taking the median of the ten highest values. The 
prominant spectral features have been identifled, following 
[Wallace et al] l|2000l ). 

The first thing to note is that the dominant spectral fea- 
tures at this resolution are due to metallic lines of Fel, Mgl, 



Si I and Til. The fluctuations in continuum level are not due 
to noise, these are blends of molecular lines. However, in this 
wavelength range and at this resolution, the lines are very 
weak and behave as a 'pseudo'-continuum. The second thing 
to note is that there is very little variation between spectral 
types; from KO to M4 the lines appear to have very similar 
strengths. This may be due in part to the narrow tempera- 
ture range spanned by RSGs (~3500-4100K), and the fact 
that all metal atoms are in the lowest ionization state at 
these temperatures. 



2.2 Model atmospheres 

For this study wc use the existing model grids of the MARCS 
project (Gustafsson ct al. 200^), which are computed in LTE 
with spherical symmetry for a range of stellar parameters. 
The mass of the star is kept as a free parameter in comput- 
ing the hydrostatic structure of the spherically-symmetric 
atmosphere, and in this grid only masses up to 5Mq are 
computed. As RSGs are typically thought to have masses 
of 15-2OM0, our use of the 5Mq masses must be justifled. 
The effect of the stellar mass in the MARCS model atmo- 
spheres is to control the so-called 'extension' (or geometrical 
thickness) ofthe atmosphere, z (for an explanation of this 
quantity see iGustafsson et al. 2008 ). Since z oc L^'^ /M, and 
2 ~ 20 in Solar units for both Red Giants and RSGs, one 
does not expect there to be much difference between the 
model spectra of stars with masses 2-2OM0 if all other pa- 
rameters are the same. Indeed, we find very little difference 
between the spectra of models with 2Mq and SM©. 

One caveat with the use of MARCS models for RSGs is 
that the atmospheres of RSGs are known to be inhomoge- 
neous and out of hydrostatic equilibrium. Convective cells 
with sizes comparable to the stellar radii and with radial 
velocities greater than the local sound speed are thought 
to permeate the star. This can lead to a patchy tempera- 
ture structure at the surface and atmospheres that are ex- 
tended (i.e. greater pressure scale heights) co mpared to that 
predicted by hydrostatic equilibrium (see e. g. Schwarzschildl 



1 19751 : IJosselin fc Plezll2007l : IChiavassa et aLlbOOgP . However, 
as the scale-height is inversly proportional to gravity, we may 
expect a model with increased scale-height to be the equiva- 
lent of a model with lower gravity. Also, as we will show later 
the metallic lines in the spectra of RSGs are very insensi- 
tive to temperature. Therefore, we expect that the effects of 
patchy surface temperature and increased scale-heights will 
have negligible impact on the derived abundance levels. 

We also note that the spectra in the MARCS model grid 
are not strictly spectra, they are samplings of the SED ffux 
at a grid of wavelengths. As the ffux between the sampling 
points is not known, if any features are 'missed' this can in- 
troduce artifacts into the spectrum when degraded to lower 
resolution (for a discussion of this effect see Plez 2008 ). We 
have run detailed comparisons between the MARCS flux- 
sampled 'spectra' and the fully-computed high-resolution 
synthetic spectrum from MARCS model atmospheres ob- 
tained from the POLLUX databasqj, and we flnd that differ- 
ences are minimal when the spectra are degraded to lower 
resolution; the only noticeable systematic difference is at the 

^ |http : //pollux ■ graal ■ uiilv-montp2 ■ f r | 
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Figure 2. An example of the fitting process for CD-31 49. Ttie 
figure shows contours of equal x^ in the log Z-T^g plane, at values 
of 0.5, 1.0, 2.0 and 4.0 above the minimum x^- The cross indicates 
the location of the x^ minimum. The gap in the contours is due 
to a missing model in the grid at Tcff=3200K, logZ=+0.5. 



blend of Fel and Si I at 1.198/ini. In future studies we will 
use fully computed spectra as they become available, but for 
now we proceed with the existing grid. 

In the MARCS grid effective temperatures are sampled 
at Teff/K = [3200,3400,3600,3800,4000,5000], which covers 
the entire range spanned by RSGs. Metallicities are sam- 
pled at log.Z=-1.5, and between -1.0 and -1-1.0 in steps 
of 0.25dex (normalized to Solar metallicity). This is easily 
sufficient to cover all metallicities expected in the Galaxy. 
The sampling in microturbulence is poorer, only values of 
^/km s~^ = [2,5] are computed. As a crude way of improving 
this sampling, we have linearly interpolated the spectra be- 
tween these models to estimate the spectra of models with 
^=3 km s~^ and 4km s~^. The sampling in log(;-space is 
also poor for these models. For now we use only the models 
with logg=(0.0, -1-1.0) (cgs units), w hich have been show n 
to be typical values for RSGs (see e.g. iDavies et al.ll2009bf ). 

The MARCS models are also available with different 
CNO abundance mixtures. The grid sampling is complete 
for the unaltered case (i.e. relative abundances of CNO set 
to Solar mixture). However, RSGs are known to display the 
results of CNO processing at their s urfaces (i.e. enhanced N, 
depleted C, see lDavies et al.ll2009bl ). For this reason we also 
examine the MARCS models with CN processing, though 
the grid sampling is poorer for these models. 



2.3 Fitting procedure 

Before comparing the observed and model spectra, we first 
degrade the spectral resolution of the model data to that of 
the observations using a gaussian convolution kernel. Then, 
in order to find the best-fit spectrum, we have experimented 
with several strategies. Firstly, we compare the normalized 



fiux at each spectral channel and compute the reduced x^ 
for each point in the model grid. 



2 

X = 






Ji, ohs Ji.i 



(1) 



where /i,obs and /i,mod are the observed and model fiuxes at 
each wavelength channel i; at is the uncertainty in the flux, 
taken to be 0.01 at each channel (i.e. SNR=100); and u is 
the number of degrees of freedom. 



N -n- 



(2) 



where N is the number of wavelength channels, and n is 
the number of model parameters being fitted. In the case of 
the MARCS grid, n = 3 {logZ, TcS, and f). We then take 
the best fitting model to be that point in the grid with the 
lowest value of x^- 

A second fitting method applies a similar principle, but 
tfiis time instead of fitting the entire spectral region we fit 
only the equivalent widths (EWs) of the strongest spec- 
tral lines: those of Fel (1.16108/im, 1.16414^m, 1.18866^m, 
1.19763/im) and Mgl (1.18314/xm, 1.20869/im). We again 
find the best fitting model by determining the point in the 
model grid with the lowest x'^- We estimate the typical un- 
certainty in EW measurements (tJi) to be around 10%, from 
repeating EW measurements of lines and slightly varying 
the wavelength ranges. 

The first 'full-spectrum' method has the advantage that 
it fits the whole spectrum, and so in principle can deal 
with blends of lines. The second 'strong-lines' method on 
the other hand takes into account only the most prominent 
spectral features and so is not swayed by unidentified blends 
for which the physical data may be poorly known. In prac- 
tice, we settled on a method which is a hybrid of these two 
- we look at the residuals between the model and observed 
spectrum, as in the 'full-spectrum' method, but only at se- 
lect regions of the spectrum. This allows us to place more 
emphasis on the strong lines, but also deal with blended 
lines such as the Fel-Sil blend at 1.198/im. 

To determine the best-fit physical properties and their 
respective uncertainties, we first find the point in the model 
grid with the lowest x^^ X'^beat- We then take into account 
all points in the grid with x'^-X'^bcst < A^^, where Ax^ is an 
arbitrary number. If there are fewer than 10 models within 
this range, we take the 10 models with the best x^- We then 
compute the weighted mean of a physical property X and 
its uncertainty aX from the relations. 



EJV 
. WiXi 
^ = ^^ — ' ^^^ 



N' j:iM^^~xr 



N' 



EN 



(3) 



where Xi is that physical property (e.g. Tofi, (,, log Z) in 
model i, and the weights Wi are calculated from the x^ oi 
each model in the grid, 



-xV2 



(4) 



and A'^' is the number of non-zero weights. The parameter 
which governs how many models are used in the computa- 
tion of the physical parameters, A^^, is initially set to 2 (see 
also Mottram et al., submitted, for similar use of this anal- 
ysis technique). However, we found that this significantly 
overestimated the errors, particularly in log Z (as will be 
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Figure 3. Examples of model fits (red dotted lines) to the data (black solid lines). Data have been resampled onto a 2x finer grid to 
help aid the identification of discrete spectral features. The top spectrum is an example of a poor fit, with reduced x^=8.9. The middle 
spectrum is a "fair" fit, with x'^=4.5. The bottom spectrum has a good fit, with x^=l-7. 



discussed in Sect. [3]), and we concluded that Ax =1 was 
more appropriate. 

An example of this process is shown in Fig. [2j where we 
show x^ contours in the Tcg-log Z plane. The cross indicates 
the derived temperature and metallicity for one star in the 
sample. It can be seen from the figure that the inferred log Z 
has only a weak sensitivity to the fitted Toff. 

As an additional check on the significance of each abun- 
dance measurement, we look at the logZ of the models in 
the grid with rcff±200K and ^±lkm s~^. In each extreme, 
we determine the 'best-fit' logZ from the model with the 
lowest x^- Therefore, we derive a range of logZ values by 
varying Tcff and ^. In the vast majority of cases, we find that 
the range in log Z is less than the sampling of the model grid 
(0.25dex), and in the worst case two grid points (±0.25dex). 
This is smaller than the typical uncertainties on each mea- 
surement (see next Section). 



3 RESULTS 

Before moving on to the derived stellar parameters for the 
programme stars, we first give examples of the quality of fits 
we achieve. In Fig. Owe show the range of fits obtained, from 
the worst to the best. In the case of the worst fit (HD 339034, 
top spectrum in Fig.[3|, it is clear that in the most part the 
strongest lines, those of Pel and Mgl, are still matched very 
well, apart from at the blue end of the displayed spectrum 
where there may be a problem with the continuum normal- 



ization. The deficiencies come mainly from unresolved and 
unidentified blends, such as that at 1.194/im, and from the 
Si I lines. However, at higher spectral resolution the Si I lines 
would be unblended and could be incorporated into the fit- 
ting procedure, possibly obtaining much better fits (as will 
be shown in Sect. 14. 2p . A thorough examination of this re- 
quires higher resolution data and is left for a future study. 
The other stars in Fig. [3] are fit very well by the proce- 
dure. This quality of fit is more typical of the results, as can 
be seen from the last column of Table [l] 



3.1 Derived stellar parameters 

We now address the results of the fitting procedure across 
all stars studied. The temperatures, microturbulent veloci- 
ties and average abundances for each star are listed in Table 
[T] In Fig. [3] we show the values of Toff and log Z derived for 
each star, as well as the averages of each parameter, as a 
function of the stars ' literature temperatu res, assuming the 
temperature scale of lLevesgue et al.l (|2005l ). In terms of tem- 
perature, we find generally acceptable agreement with the 
literature values. We do find on average that the tempera- 
tures we derive are cooler by «150K, with a la standard de- 
viation of 170K on the distribution. This effect is more pro- 
nounced at higher tempera tures. However, if one takes th e 
old temperature scale from iHumphrevs fc McElrovi ([198J) , 
which are on average ~200K cooler than those of Levesque 
et al., the agreement is almost perfect. This is curious, as 
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T able 1. Derived param eters for the programme stars. Temperatures are those appropriate for their spectral types from the calibration 
of lLevesgue et alj l|200 5fl: where a star is between spectral types we linearly interpolate between the two. Abundance types are either 'st' 
for standard, or 'CN' for CN processed. 



Name 



Spec Type T(K) Tfit(K) 



log[Z] 



^/km s 



AbT 



HD63302 


Klla-Iab 


4100 


3660 ± 180 


0.02 ± 0.20 


3.4 ± 0.5 


St 


4.47 


HD212466 


K20-Ia 


4015 


3770 ± 170 


0.17 ± 0.20 


3.8 ± 0.5 


st 


5.57 


HD187238 


K3Iab-Ib 


4015 


3630 ± 180 


0.04 ± 0.21 


2.3 ± 0.5 


st 


4.91 


HD185622A 


K4Ib 


3928 


3590 ± 170 


0.12 ± 0.20 


2.4 ± 0.5 


st 


5.95 


HD216946 


K5Ib 


3840 


3580 ± 230 


-0.15 ± 0.27 


2.5 ± 0.5 


CN 


3.29 


HD14404 


Ml-Iab-Ib 


3745 


3580 ± 170 


0.10 ± 0.19 


2.5 ± 0.5 


CN 


4.45 


HD39801 


Ml-M2Ia-Ia 


3710 


3520 ± 160 


0.19 ± 0.21 


2.3 ± 0.5 


st 


6.59 


HD35601 


M1.5Iab-Ib 


3710 


3630 ± 160 


0.11 ± 0.19 


2.8 ± 0.5 


st 


3.87 


BD+60 265 


M1.5Ib 


3710 


3570 ± 140 


0.16 ± 0.17 


2.3 ± 0.5 


st 


6.90 


HD339034 


Mlla 


3745 


3820 ± 160 


0.29 ± 0.15 


3.5 ± 0.5 


st 


8.86 


HD 10465 


M2Ib 


3660 


3520 ± 220 


-0.08 ± 0.27 


3.4 ± 0.6 


st 


1.78 


HD14469 


M3-M4Iab 


3550 


3530 ± 160 


0.15 ± 0.22 


2.4 ± 0.5 


CN 


5.76 


CD-31 49 


M3Iab-Ia 


3605 


3580 ± 220 


-0.15 ± 0.29 


3.1 ± 0.6 


st 


1.67 


RWCyg 


M3toM4Ia-I 


3550 


3680 ± 180 


0.31 ± 0.21 


2.6 ± 0.5 


st 


6.98 
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Figure 4. Results of the analysis process on the programme 
stars, as a function of their temperatures inferred from their lit- 
erature spectra l type s and the temperature calibration scale of 
iLevesaue et al.l 1(20051) . The top panel shows the difference in Toff 
between the literature values and that derived by the fitting pro- 
cedure. The bottom panel shows the derived average metal con- 
tent log Z. The weighted means of each parameter arc shown in 
the panels. 



the RSG temperature scale was rederived from the current 
generation of MARCS model atmospheres, but concentrat- 
ing mainly on the R and I spectral windows. We will discuss 
how the derived temperature may be affected by line blend- 
ing in Sect, [l] but we will also show in Sect. 13.21 that the 



fitted stellar temperature has very little impact on the mea- 
sured abundances, so for now we are not concerned by this 
minor discrepancy. 

The fitted values of log Z are very tightly distributed 
around zero, with mean metallicity of logZ=0.13 ± 0.14, 
where the uncertainty again is the la standard deviation. 
These values are therefore consistent with Solar metal abun- 
dances, which is a reassuring result for a sample of objects in 
the Solar neighbourhood. If we assume for the moment that 
all objects in the sample have log Z=0.0, then the bottom 
panel of Fig. |4] serves to illustrate an empirical determina- 
tion of the precision of our analysis method, that is that 
we are deriving abundances accurate to ±0.14dex. However, 
there is likely some intrinsic spread in the metallicities of 
the objects in our sample, so this value is an upper limit to 
the experimental uncertainty. This compares to the average 
errors we determined in the fitting process of ~0.2dex (see 
column 4 of Table [l]). So, it is likely that our analysis pro- 
cess overestimates the experimental errors slightly. When 
we used Ax^~2 in the fitting process (see previous Section) 
these errors were higher still, at around ±0.25-0.3dex, yet 
the average log Z and its standard deviation were the same 
as with Ax^ = l. This is the reason we opted for Ax^~l 
during the fitting, since the errors were more comparable to 
'empirical' error of ±0.14dex. This level of precision is as 
good as any other abundance measurement technique. 

For the microturbulent velocities ^, all measured values 
are between 2-4 km s~^, with an averag e of 2.8 ±0.5 km s~^ . 
This is a typical va lue for RSGs (see e.g. lPavies et al.ll2009bl : 
ICunha et al.ll2007^ . We note however that the ^-sampling of 
the MARCS grid is very coarse, with computed models at 
2 and 5 km s~^ only. Our grid points at 3 and 4 km s~^ 
were simple linear interpolations in between the computed 
models, and so if non-linear effects are present there will be a 
systematic uncertainty in our derived values of ^. However, 
we can say that the current ^-sampling of the model grid 
encompasses the range of values expected from the RSGs in 
our sample. The impact of the uncertainty in ^ on abundance 
measurements will be explored in the next Section. 

We explored two sets of models with different CNO mix- 
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Figure 5. The impact of fixing various input parameters on tlic 
results of the fitting process. Each panel shows the results of the 
best fit (i.e. letting all parameters vary), fixing ^=2 km s~^ and 
^=4 km s~^, fixing the CNO elemental abundances to the 'stan- 
dard' and 'processed' values, and fixing logg=-|-1.0. The top panel 
shows the weighted-mean abundance level of all stars and its stan- 
dard deviation. The bottom panel shows the median x^ of all fits 
from each of the fitting methods. 



tures - the 'standard' set, and the CNO-processed set with 
depleted C and enhanced N. For the majority of the objects 
studied the standard CNO abundance models provided the 
best fits. As with ^, our sensitivity to this free parameter 
will be explored next. 



3.2 Robustness of abundance measurements 

In the previous Section we have shown that the abundances 
we measure in this technique are extremely consistent with 
the expected values, with a very high level of precision. To 
investigate the robustness of these measurements, we now 
explore the effect of varying the free parameters on the re- 
sults. In Fig. Owe plot the weighted mean of all stars' fitted 
abundance levels, the standard deviation on this value, and 
the reduced x^ values when forcing each of the model input 
parameters to be fixed. 

Sensitivity to microturbulence This parameter affects 
the EWs of strong lines. If lines are saturated, then by 
adding a degree of microturbulence the overall capacity to 
absorb is increased, and stronger absorption lines are ob- 
served. The second and third points on the y-axis in Fig. 
[S] show the impact of forcing the fitting routine to find the 
best fit when the microturbulence is fixed at ^=2 km s~^ and 
^=4 km s^^. A systematic effect on the derived \ogZ can be 
seen, that is for lower ^, higher log Z is found. However, 
the difference is small (±0.2dex), while larger dispersion in 
log Z is found, illustrated by the size of the error bar. The 



larger dispersion indicates that, on average, poorer fits are 
found, in accordance with the larger values of x^ from the 
fitting procedure (bottom panel of Fig. [5} . The small change 
in the average log Z suggests that the uncertainties in abun- 
dance levels propagated through from uncertainties in ^ are 
typically less than the experimental error in log Z. 



Sensitivity to CNO abundances We find very little sen- 
sitivity to this input parameter. Whether the 'standard' or 
'CNO processed' abundances are used, the derived log Z is 
very close to the best-fit value, with similar error-bar (see 
the fourth and fifth points on the y-axis in Fig. [S}. 



Sensitivity to gravity It is not possible as yet to make 
as thorough an investigation of the effect of varying this 
parameter as for the others, as it is poorly sampled in the 
MARCS grid for supergiants. However, models do exist for 
\ogg=+1.0 and 0.0. These values probably bracket the ac- 
cepted values for most RS Gs, with logg=0.0 being the most 



commonly used number (ILambert et al. 1984 : ICarr et al.l 
I2OOOI : iLevesaue et"al]|2005l : ICunha et al.ll2007l ). Though val- 
ues of logg as low as -0.5 have been measured for some 
RSGs, the vast majorit y of RSGs have gra vities of loggr^O.O 
and higher (see Fig. 5 of lCunha et al.|[2()07l ). We find that, on 
average the logg^-l-l.O fits produce significantly larger x^, 
and a slightly larger spread on the log Z distribution. The 
average log Z however remains consistent with the best-fit 
solution to within the errors. We experimented with using 
the log(7=-|-1.0 models in the formal fitting procedure, and 
deriving best-fit values for log g rather than assuming a blan- 
ket value of log g=0.0. In all cases we found that the derived 
values of logg were between 0.0 and 0.2, with negligible ef- 
fect on the derived logZ and Teff. 



Sensitivity to input Toff The top panel of Fig.[5]illustrates 
that, for any set of model parameters used, we consistently 
derive temperatures that are systematically lower than the 
literature values, especially for earlier spectral types. As 
commented upon earlier, this may be due in part to the 
recently-rederived temperature scale for RSGs. Neverthe- 
less, to investigate the significance of this effect we re-ran the 
fitting procedures but this time forced the fitted tempera- 
ture to be the same as the literature temperatures, allowing 
the rest of the parameters to vary. The results are shown in 
the right-most points of Fig. \E\ Though the reduced x^ val- 
ues are slightly worse than for any of the other model runs, 
the average abundances are again extremely consistent with 
the best-fit model, only with a slightly larger dispersion. 

From this analysis, it appears that the free parameters 
the fitting process, namely T^g, ^, logg, and X(CNO), have 
very little effect on the derived abundances. The standard 
deviation on the results of each model run in Fig. [5] is only 
0.2dex, with a mean of -|-0.0dex. These numbers are com- 
parable to those derived for the 'best-fit' results. Thus, it 
seems that any systematic uncertainty in the fitting process 
is only comparable to the random error, and the derived 
abundance level is robust to fluctuations in the model's free 
parameters. 
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Figure 6. The effect of average metal content log Z on the EWs of 
the diagnostic lines used in this study, as measured from the syn- 
thetic spectra in the MARCS model grid. Points show the average 
EWs, while error bars show the range of EWs measured across 
the effective temperature range of RSGs (3400-4000K). All mod- 
els use log(;=0.0, 5=3km s""*^, and standard CNO abundances. 
The X- values of each of the curves have been slightly offset from 
one another for clarity. 



4 DISCUSSION 

4.1 RSGs as abundance probes 

The results of this study suggest that the MARCS model 
atmospheres provide excellent fits to the J-band spectra 
of Solar-neighbourhood RSGs. Moreover, the derived abun- 
dances seem to be relatively insensitive to the model input 
parameters Teff, logg, ^ and the mix of CNO abundances. 
Therefore, by constraining these input parameters to the 
ranges that are known to be valid for RSGs, it should be 
possible to derive metal abundances for RSGs to a precision 
of O.ldex using medium resolution J-band spectroscopy. 

To illustrate the sensitivity of the J-band diagnostic 
lines to average metal content, in Fig. [5] we plot the strengths 
of the strong lines identified in Sect. l2.3l as a function of log Z 
for the model spectra in the MARCS grid. The error-bars in- 
dicate the minimum and maximum EWs observed across the 
range of RSG temperatures, from 3400 to 4000K. All mea- 
surements were computed from the \ogg—Q.Q, ^=3 km s~^ 
models with standard CNO abundances. The plot shows 
that, for these lines, there is very little variation with 
RSG spectral type, presumably due to the narrow range 
of temperatures that RSGs span. By contrast, the effect 
of changing metal content is much more pronounced. For 
the typical abu ndance gradient in a spiral galaxy of 0.05- 
0.2dex/kpc (s ec Garnctt et al. ' 19971 : iKudritzki et al.ll2008l : 
iBresolin et"lj.i 2009; U et al. 200i), the strengths of these 
hues could increase by ~50% between the innermost and 
outermost regions. As such variations in line strength would 
be easily detectable, RSGs could be powerful tools with 
which to study abundance trends in their host galaxies. 



4.2 Benefits of increased spectral resolution 

The spectra studied here were taken from the IRTF spec- 
tral library, the spectral resolution of which is R = 2000. 
At this resolution we have identified six spectral lines which 
have sufficient strength to be largely unaffected by blending. 
Nevertheless, there is clearly some blending present between 
the metallic lines, e.g. those of Fel and Til. At lower metal- 
licities, where the metallic lines will tend to have lower EWs, 
blending may begin to introduce problems. 

However, at greater spectral resolution the blending of 
the major lines is greatly reduced, increasing the number of 
strong diagnostic lines that may be studied, as well as in- 
creasing the number of elements that can be probed. This is 
illustrated in Fig. [T] where w e plot the spect r a of a Her and 
a Ori originally presented in lWallace et al.l (|2000)- At this 
spectral resolution of _R = 6000 we see the lines of Si I and 
Til as unblended from the Eel lines, particularly at 1.189/^m 
and 1.198/im. The results of the fitting process for these two 
stars are shown in Table [21 Note that the abundance mea- 
sured for Q Ori of log Z =+0.24^: . 22 is entirely consistent 
with the results of iLambert et all (|l984 Fe/Hsi+O.l) 

The benefits of this are twofold. Firstly, it is possible to 
add the Si I and Til lines to the list of 'strong features' that 
are measured during the fitting process, enabling a more ro- 
bust fit and determination of stellar temperature. For exam- 
ple, using only the Fe I and Mg I lines in the fitting procedure, 
we measure a temperature for a Ori of Teff=3300±200K. 
However, these fits underestimate the strengths of the Si I 
lines, as in the case of HD 212466 (see Sect. |3}. When the 
Si I and Til lines are folded into the analysis we measure 
Jcft=3660±170K, which is much closer to the temperature 
appropriate for its spec tral type (Ml. 5 la b, rcff=3710K, us- 
ing the calibration of Levesque et al.ll2005l ). We note however 
that the derived metallicity is relatively insensitive to this 
effect, changing by only ~0.1dex. 

Secondly, at increased spectral resolution the number of 
metallic elements that it is possible to investigate is greater. 
By studying the combined abundances of Mg, Si and Ti, 
combined with more detailed modelling, it should be possi- 
ble to make a quantitative comparison of the relative abun- 
dances of Fe-group and a-group elements, a key diagnostic 
of a galaxy's star-forming history. Though the spectra we 
show in Fig. [7] have resolution of _R = 6000, we have found 
that it is still possible to separate the spectral lines at reso- 
lution of only R = 3000. As lower resolutions mean shorter 
exposure times for objects of a given brightness, the limit- 
ing distance at which RSGs could be studied in this way is 
increased. 



4.3 Extragalactic potential 

The potential of RSGs as extragalactic abundance probes is 
enormous. With absolute magnitudes of Mj ^ —8 to —11 
mag and very characteristic colors ( 1/ — 7 ~ 2.0 mag) they 
can easily be detected as individual objects in galaxies at 
large distances. Using new IR multi-object spectrographs 
such as KMOS at the ESO VLT or MOSFIRE at Keck one 
will be able to reach limiting magnitudes of Mj fa 19 mag 
with a spectral resolution of _R=3000 and a signal-to-noise 
of 100 in long exposure sequences (for instance, totalling 
two nights). In this way quantitative spectroscopy of RSGs 
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Table 2. Derived parameters for the two stars for which higher spectral resolution was available. 
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Figure 7. Best fits to the J-band spectrum of a Her and a Ori, with spectral resolution R = 6000, published in lWallace et al.l 1120001 ). 
See text for details. 



out to 10 Mpc distances will be possible allowing precise 
chemical abundances of galaxies to be determined and the 
evolution of galaxies to be disentangled. 

The future, with next generation of 30m+ ground-based 
telescopes and the James Webb Space Telescope, is even 
more promising. These new facilities will be optimized for 
very high spatial resolution observations in the IR, so that 
the gain in limiting fluxes for multi-object spectrographs 
increases with the fourth power of aperture diameter. For 
instance, using the E-ELT exposure time estimator and in- 
cluding the dramatic effects of adaptive optics we will be 
able to go down Mj ~ 22 mag, which will give us a limiting 
distance of 30 Mpc. This opens up a truly remarkable vol- 
ume of the local universe for very detailed and precise stud- 
ies of the formation and chemical evolution of entire galaxy 
clusters, such as Virgo, Fornax, Puppis and Eridanus. 



other abundance indicator. We find that the derived abun- 
dance levels are relatively insensitive to the other model pa- 
rameters of effective temperature, microturbulent velocity, 
CNO mixture and surface gravity, though the effects of very 
low surface gravities (logg<0.0) have not been investigated 
here. 

Though we do find effective temperatures which are on 
average lower than the literature values by ~100-200K, this 
has very little impact on the derived abundances since log Z 
is extremely insensitive to Toff. Forcing the fitted tempera- 
ture to be the same as the literature temperature for each 
star resulted in output abundances that were very similar 
and within the errors of the best fits. Furthermore, we ar- 
gue that this Tcft discrepancy may be solved at resolutions 
of 7? > 3000 when the blending of key diagnostic lines is 
diminished. 



5 CONCLUSIONS 

To demonstrate that it is possible to get accurate and pre- 
cise abundance measurements from Red Supergiants (RSGs) 
at modest spectral resolutions of _R= 2-3000, we have stud- 
ied a sample of RSGs in the Solar neighbourhood in con- 
junction with the latest MARCS model atmospheres. With 
no a priori assumptions about stellar parameters, we derive 
chemical abundances for each object which are close to the 
Solar value, with a mean abundance for the whole sample 
of 0.13±0.14dex. This dispersion in the mean is notwith- 
standing any intrinsic spread inherent in the sample, and is 
indicative of our experimental uncertainty. The level of pre- 
cision achieved by this method is therefore as good as any 



The results of our work therefore suggest that it is pos- 
sible to obtain accurate abundance information at spectral 
resolutions which are much lower than commonly used at 
present. This opens the door for abundance studies of entire 
galaxies at infrared wavelengths, out to distances of several 
Mpc. 

As for the future, we will shortly present an analysis of 
a sample of RSGs from the two Magellanic Clouds, showing 
that this technique can be applied across a large range of 
metallicities. We are also in the process of expanding the 
model grid to include a broader range of gravities, better 
sampling in ^-space, and with accurate synthetic spectra 
computed at high resolution. 
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